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A model for autonomous feedback control of particle transport through a large number of channels 
is introduced. Interactions among the particles can lead to a strong suppression of fluctuations in 
the particle number statistics. Within a mean-field type limit, the collective control mechanism 
becomes equivalent to a synchronization with an external clock. The diffusive spreading of the 
feedback signal across the channels shows scaling, can be quantified via the flow of information, and 
shows up, e.g., in the spectral function of the particle noise. 
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I. INTRODUCTION 

Feedback loops are an interesting tool to modify and 
control the transport dynamics of systems, both classical 
and quantum [iHin]. Increasing efforts have been made 
recently to understand and quantify closed loop control 
schemes (e.g. in Maxwell demon, information tape or 
network models) from the perspective of thermodynam¬ 
ics, statistical mechanics and system-bath theories ffn- 

127]. 

A key question here is how to design the feedback loop 
itself, i.e. whether to regard it as part of some active ex¬ 
ternal measurement scheme, or ‘passively’ as an extended 
part of the system itself. This latter form of coherent con¬ 
trol [751 - I55| is particularly appealing for quantum systems 
(also cf. [35] for coherent feedback control of quantum 
transport and further references), as it avoids the need 
to involve the measurement process. But also classically, 
a full microscopic understanding of feedback control re¬ 
quires the modeling of some form of interaction between 
the system and its controller. 

It is this perspective from which we re-visit a feed¬ 
back control scheme [32] of a stochastic process describ¬ 
ing (source to drain) transport of particles, which in its 
simplest version describes a continuous, active modula¬ 
tion of transition rates conditioned upon the full counting 
statistics (FCS) |3S], i.e. the statistics p(n, t) of particle 
numbers n that have accumulated in a drain reservoir 
after a certain time t. In the new feedback model, this 
modulation is achieved in passive mode via interactions 
among particles in N coupled transport ‘channels’ that 
provide the feedback by making the individual transi¬ 
tion rates in each channel dependent upon the (reservoir) 
state of the entire system. 

The starting point in this paper thus is an infinite sys¬ 
tem of coupled Poissonian processes for the FCS obtained 
by adiabatic elimination of internal ‘connector’ degrees 
of freedom between source and drain reservoirs, cf. Fig. 
(la). Detailed balance for the transition rates then al¬ 
ready restricts the form of possible feedback models that 
can be obtained from microscopic interactions. Choosing 
a classical limit of weakly interacting uncharged fermions, 
we interpret the interaction-induced strong reduction of 
fluctuations as a ‘condensation’ in the space of reservoir 


particle numbers. This is accompagnied with an instan¬ 
taneously increased flow of information between a given 
channel and the rest of the system. For periodic channel 
structures, Fourier analysis yields detailed predictions for 
the diffusion-like spreading of the autonomous feedback 
signal across the channels, which can be made visible in 
observables like the particle noise spectrum. 


II. MODEL 

Our model is defined by a d-dimensional lattice with 
N lattice sites R;, I — 1, each of which serves as 

a connector within a channel with particle flow from a 
source to a drain reservoir at rate 7;+ (‘forward’), and 
backwards from drain to source at rate 7;_, cf. Fig. (la). 
There is no direct transfer of particles between different 
channels, but the time-independent rates ■ji± = 7/±(n) 
are allowed to depend on the state n = (rii, ...,n7v)^ of 
all the drain reservoirs. These are defined by the numbers 
n; G Z of additional particles in drain reservoir I, when 
counting starts at time t = 0. 

The dynamics is described by a Markovian rate equa¬ 
tion for the probability p(n, t) of the drain reservoirs to 
be in state n at time t, 

P(n, n')p(n', t) - r(n', n)p(n, t)], (1) 

n' 

where r(n',n) = 7i±(»)<)„',„±ei with the 

Cartesian unit vector e/ = (0,..., I, ...0)^, such that we 
can write 

N 

1 = 1 ,± 

which makes the contributions from forward and back¬ 
ward jumps in the N channels explicit. As initial condi¬ 
tion for p{n, t) we set p(n, 0) = (5„^o- 

The rate equation Eq. ([^ can be generalized to the 
more complex situation where the N connectors between 
the source and drain reservoirs have themselves an in¬ 
ternal structure, i.e. internal states (like in quantum 
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FIG. 1: (a) Particle transport between source and drain 

reservoirs at rates 7 ± through connectors in channels at po¬ 
sitions Ri. (b) Connectors could have internal states, (c) 
Autonomous feedback via interactions (g > 0) among drain 
reservoirs reduces fluctuations of particle numbers ni, repre¬ 
sented as Brownian particles on a line. This is compared to 
the non-feedback case g ~ 0 during time evolution from to to 
ti. (d) Active, non-autonomous measurement-based feedback 
for a single channel |34| . 


For example, a linear modulation f±{x) = 7±(1 ± gx) 
with control parameter 0 < (; <C 1 continuously compen¬ 
sates particle number fluctuations in that the forward 
rate 7 + is decreased if there are too many particles as 
compared to a reference charge It. The rate equation 
can then be solved exactly, and the main result [34] is 
a freezing of the single channel full counting statistics 
Pi{n,t) at large times t around a mean (n)* oc t. Quan¬ 
titatively, the second cumulant follows from Eq. ([^ as 

= (4) 

which means that fluctuations are strongly suppressed 
for g > 0, which is schematically shown in Fig. (fc). In 
a similar way, expressions for higher cumulants [34j can 
be derived that all saturate at finite values oc 1 /g at large 
times t. 


B. Interacting channels 

Alternatively, feedback occurs via the design of inter¬ 
actions among particles. In most of what follows we will 
be dealing with an autonomous form 


dots) that can be occupied by particles. It seems diffi¬ 
cult to make progress then (apart from large-scale numer¬ 
ics) , but on the other hand the above rate equation then 
can follow as an effective model via adiabatic elimination 
when time-scale separation is possible. For example, this 
is the case with one internal state connected to a source 
at a very fast rate and a drain at a very slow rate (not 
shown here). In all what follows, we will therefore work 
with Eq. ([^ as the starting point. 

The rates 7 i±(n) are the key elements in describing 
the interactions among the channels and the effective 
feedback mechanism created thereby, which is why we 
discuss them in quite some detail in the following. 


A. Decoupled channels with time-dependent rates 

In the simplest version of the feedback model, the rates 
have the time (t) and n-dependent form 

7 i±(n) = /±(-f;i-«-/), (3) 

where the function f±{x) describes a synchronization be¬ 
tween the actual number ni with an external reference lit 
that grows linearly with time (‘clock’ with a fixed control 
current Ii). The form Eq. (|^ then immediately reduces 
Eq. <§ to decoupled rate equ ation s for the single-channel 
probabilities Pi (n, t), cf. Eq. ( |A4[ ) in Appendix A. Physi¬ 
cally, this kind of feedback can be achieved via an instan¬ 
taneous control, i.e. a modification of the rates ‘by hand’ 
in an active way by an external agent like an electronic 
circuit, cf. Fig. Id). 


li± (n) = 7± 


N 

-ni) , 
i=i 


(5) 


which is parametrized by interaction matrix elements gij 
and which is linear in the ni. The rates Eq. ([^ are a 
generalization of the Eq. ([^ in the following way: con¬ 
sider the simplest all-to-all coupling gij = for which 
the rates have the form 7 ;± (n) = f± — ni), with the 
average riav = (ni + ... + n]sf)/N and f±ix) = 7±(1 ± 52 ;) 
as above. Now, for A^ —)■ 00 the average riav becomes 
a macroscopic variable, and neglecting the fluctuations 
of the latter against those of n; one obtains the rates 
Eq. ([^ with lit = (riav), he. we recover the synchroniza¬ 
tion model in this mean-field limit (see Appendix A for 
some more details). 

The rates Eq. have a very intuitive form: the n/- 
dependence leads to a dynamical compensation of parti¬ 
cle number differences among the channels. If there are 
too many particles in the drain reservoir I (as compared 
to all the other drains j ^ 1), the forward (backwards) 
rate for channel I is decreased (increased), and the other 
way round if there are too many particles in drain 1. As 
we will show below, this compensation mechanism tends 
to suppress fluctuations of the drain reservoir particle 
numbers ni also within each single channel 1. 

The rates Eq. ([^ can potentially become negative (and 
the model unphysical) in the linear approximation; how¬ 
ever, this is not really a problem as long as the interac¬ 
tions are sufficiently small. More importantly, with the 
gij as unspecified parameters, we have not yet addressed 
the issue of detailed balance between forward and back¬ 
ward rates so far. 
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C. Detailed balance 

We therefore consider a dependence of 7 i+(n) = 
ri(£i(n)) on the energy differences ei(n) = A/r; + A£’;(n), 
where Afii = is the difference of the chemical 

drain and source potentials in channel I, and AEi{n) = 
E{n + ei) — E{n) with the energy E{n) due to interac¬ 
tions among the particles in all drain channels. Corre¬ 
spondingly, 7 i-(n -I- ei) = ri(—£i(n)), and detailed bal¬ 
ance reads 


7i+(n) = 7i_(n-f ei)e / 3 Aw-i-A£,(n)]^ 


( 6 ) 


where all the channels are kept at the same temperature 
= /?-!. 

By expanding the rates r/(£i(n)) and the r.h.s. of 
Eq. Q in the interaction energies A£'/(n), we find a lin¬ 
earization 

7 ,±(n) = 7 ;±x (^l + ^|||^[£;(n±eO-E(n)]) (7) 

with 7 jd- = ri(±Afii), which is consistent with detailed 
balance if AEi{n) is sufficiently small. 

We can be more specific in Eq. ([^ if we model the par¬ 
ticles as fermions that tunnel between a non-interacting 
(source) and an interacting (drain) region, similar to 
the orthodox model of Coulomb blockade for charged 
fermions (electrons) in metallic quantum dots [37j . In 
this case, the specific expression 


riis) = Ti- 


./3e _ 1 


( 8 ) 


can be derived, where ri is a parameter that depends 
on the (fiat) density of states in source and drain reser¬ 
voirs. The coefficient in the rate Eq. 0 then becomes 

n(±Aij 7 ~ ~/3/2 provided /3|A^i| ^ 1, in which case 
detailed balance is always fulfilled in lowest order in 
f5AEi{n) in Eq. ([ 6 ). This also shows that in this ex¬ 
pansion, we are dealing with rates in a high-temperature 
limit, where the particles essentially behave as classical 
objects and quantum statistics and quantum coherence 
plays no role. What we are left with as a key modeling 
parameter are the interactions among the particles. 


D. Interaction form 

In order to make analytical progress with Eq. we 
further linearize the rates 7i±(n) in n by assuming a 
quadratic dependence 


£;(n) = 


( 9 ) 


with a symmetric N x N matrix U of interaction param¬ 
eters whence Eq. 0 becomes 


7 /±(n) = 7 /± 


I T |e/ Un 


( 10 ) 


where we neglected a constant, n-independent term. 

From Eq. ( [T^ we recognize that we obtain the feedback 
form Eq. ([^ with the identification 


N 


Uu = 


Sij ^ 9ii - 


9ij 


( 11 ) 


and homogeneous rates 7 z± = 7 ± independent of channel 
index 1. For example, with gij = ^ and a single interac¬ 
tion parameter g > 0 , one has Uij = ^g{Sij — which 
mimics a mean-field type repulsive interaction within the 
same channel I = j, and an attractive interaction between 
particles in different channels I ^ j, regardless of particle 
distance. The form Eq. (Ill means that the interaction 
matrix U has a zero eigenvector (l,!,...,!)"^ such that 


N 


E^b = o 


( 12 ) 




and the (repulsive) interaction strength Uu within each 
channel I is cancelled by the sum of the remaining (attrac¬ 
tive) interactions with all other channels j. Physically, 
such an interaction corresponds to a repulsive short-range 
and an attractive long-range, van-der-Waals like force be¬ 
tween uncharged particles. 

The choice of parameters in the interaction matrix Uij 
in Eq. 0 is restricted to the particular form Eq. ( [TT] ), 
because we demand a stationary particle current to flow 
through the channels, i.e. a situation where the average 
numbers n; continue to grow linearly with time t. 

We obtain the equation of motion for these averages 
(ni) = from Eq. 0 by summation over 

as 


— (n;) = (7i+(n) -7,_(n)) 


(13) 


or |(n) = A7 + (G++G-)(n), in matrix form with the 
matrix G)^ = and the vector A-Ji = 7 /+ — 7 i_. 

For general Uij, this equation has a fixed point with (n) = 
— G“^A7 for invertible G = G+ -I- G~^. Here, however, 
we are interested in the opposite case where G (and thus 
U) is not regular but has a zero eigenvector (l,!,...,!)"^ 
leading to Eq. (12) and Eq. (II). 


As an immediate consequence, we find the solution of 


Eq. (13) 


im) = ( 7 + - 7_)t, 


(14) 


i.e. a linear increase of the particle numbers in all chan¬ 
nels regardless of the interaction strength. 


III. FLUCTUATIONS 

The average flow of particles has the rather trivial form 
Eq. (14) for the homogeneous case considered here and 
in the following, i.e. equivalence of all channels I with 
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translational invariance across the d-dimensional lattice 
of connectors between source and drain reservoirs. As in 
our previous work [33], we expect the feedback to dras¬ 
tically modify the full counting statistics of the particle 
numbers ni. The most important object to quantity fluc¬ 
tuations then is the second cumulant. Higher cumulants 
are difficult to extract from Eq. ([^, except for special 
cases as in [33] , and we will therefore effectively evaluate 
only a Gaussian model for the fluctuations. 


A. Second cumulant 

We quantify the fluctuations by the correlation func¬ 
tion 


Cu'{t) = (ninv) - {ni){nv), (15) 

for which we obtain an equation of motion via Eq. ([^. 
Multiplication of Eq. ([^ with (nei)(ne//) and summation 
over n yields ftinm,) = J2±{Su'li±{^)Tnni'±{n) + {l o 
I'), and for the homogeneous case with rates Eq. ([^ we 
obtain 


d 

dt 


Cw 


= 7 



N 


+ X! 3i'j 

i=i 


Cw) + l^l 


(16) 


with the definition 7 = 7 + -|- 7 _. We can now take advan¬ 
tage of the discrete translational invariance on the lattice 
{R/} of channel connectors I by defining the Fourier de¬ 
composition Cw = ^ X^kk' in recipro¬ 

cal space for the correlation function and 




,-jk(Ri-Rj 


^s(k) 


(17) 


for the interaction matrix elements in Eq. ([^. 

We solve Eq. (16) for an initially empty drain reservoirs 
with C'n'(O) = 0, 




(18) 


The k = 0 contribution in the sum Eq. (18 1 is extracted 
with the limit Tk^o = 0 ; and we have 






,-279ii 


( 21 ) 


where we used X^k^o ® **^^*^" 
second term. 


1 H-i') = NSiii — 1 for the 


The result Eq. (21) already captures one of the main 


ingredients of the feedback mechanism described by our 
model: at large times t, the fluctuations Cw{t) increase 
linearly but are suppressed by a factor of 1/N as com¬ 
pared to the simple non-interacting case with the Pois- 
sonian fluctuations {6nf{t)) = yt. 

In the limit of a very large number of channels, for 
times such that iV 3> yt S> 1 we obtain a freezing of the 
fluctuations towards the value Cw{t) Sw 1‘^g- For the 
fluctuations within one channel, this becomes visible as 
a plateau at intermediate times in Fig. ([^), where we 
plot Cii{t) for various values of N. 

In each channel I, a synchronization occurs where in 
the rates Eq. the contributions from all the other 
channels can be replaced by a time-dependent average 
Eq. (§ with lit = ( 7 + — 7 _)t. This is in accordance |84] 
with the non-autonomous model Eq. ([^ , where we found 
a freezing of the entire full counting statistics around a 
moving mean value Eq. (141 at large times t, with a sec¬ 
ond cumulant Eq. (|^ that co-incides with Eq. (21) for 
N —>■ 00 . In our mean-field model here, instead of syn¬ 
chronization with an external ‘clock’ current Ii, we have 
auto-synchronization between all channels. In Appendix 
A, we demonstrate the equivalence between the synchro¬ 
nization model Eq. and the A-channel model with 
infinite range interactions in the mean field limit by de¬ 
riving the equation of motion for the single channel prob- 
abibilty pi{n, t). 

Diffusion limit and scaling. — We now evaluate the re¬ 
sult Eq. (18) in the opposite limit of short-range inter¬ 


actions gij, which has a more interesting dynamics as 
compared to the mean-field case above. We can evaluate 
the Fourier-matrix elements g(k) for nearest-neighbor 
interaction in Eq. (jll 


with the definitions 




— Tki 


Tk = 27(g(0) - g(k)), (19) 


which are real quantities as we recognize from the repre¬ 
sentation g(k) = the symme¬ 

try gw = gi'i- 

Infinite range interaction . — We first evaluate the re¬ 
sult Eq. (18) for the mean-field type model gij = in 
the interactions Eq. which (simple as it is) already 
contains the essential ingredients of the feedback mecha¬ 
nism we try to analyse here. 


Using g(k) = ^ E/z'e 

Eq. (jl^ we obtain 


ik(Rz—Rj/) _ 


= g6k,o, from 


( 20 ) 


gij=gS(^ij), g>0. 


( 22 ) 


For example, on an infinite lattice with lattice constant 
a in d = 1 dimension, we have g{k) = 2g cos ka. 

The decay rates Eq. (19) acquire the form 


Tk = 275 a^k^ 


0 


(23) 


in the continuum limit where the connector lattice con¬ 
stant a approaches zero. In the terminology of dynamical 
phase transitions, this corresponds to model-A dynam¬ 
ics [381 ESI due to the fact that the total particle number 
Wot = E/ the drains grows with time and is not 
conserved. 

The diffusive nature of the feedback dynamics is now 
clear from the second cumulant Eq. (18), the temporal 


Fk^^o = Sgy. 
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FIG. 2: (a) Diagonal second cumulant Cu Eq. (211 multi¬ 

plied with interaction parameter g and (b) information cur¬ 
rent Eq. (401, as a function of time t multiplied wit h 

pq. (c) Stationary diagonal noise spectrum (Z = Eq. (25l) 
as a function of scaled frequency fZ = Lo/{gj). Curves for in¬ 
finite range interaction model Eq. ( |20[ ) with channel numbers 
= 50 (magenta dotted), N = 500 (red dashed), N = 5000 
(blue thick line), diffusive approximation Eq. (23l in d = 2 
dimensions with finite N in (a), and nearest neighbor model 
Eq. (221 in d = 1 dimension with N ^ oo (thin black line). 


derivative of which. 




„-ik(Ri 


-R,,)g-2nk2t 


(24) 


has the form of a d-dimensional diffusion propagator, 
where D = jga^ is an effective diffusion constant. Ac¬ 
cordingly, we obtain the simple, d-dependent temporal 
and spatial (with respect to |Ri — Rpl) scaling behavior 
of Cw (t) characteristic of a critical dissipative Gaussian 
model |39j with dynamical critical exponent z = 2. 


For any finite number N of channels, the integral ap¬ 
proximation to the k-sum Eq. (18 1 has a lower cutoff 
at the inverse system size oc N~‘^. With increasing d, 
at intermediate times Cii{t) develops a flat plateau as 
in the all-to-all coupling model Eq. (20 1 , before the lin¬ 
ear in t behaviour sets in at very large times. This is 
shown in Fig. ([^) for the diffusive model Eq. (231 in 
d = 2 (with the same values for N as in the all-to-all 
coupling model), and for the full nearest-neighbor inter¬ 
action model Eq. (22 1 in d = 1 in the limit N = oo, where 
Cii{t —>■ oo) oc 


B. Noise Spectrum 

The above temporal fluctuations can be related to the 
fluctuations 5Ii{t) = Ii{t) — {Ii) of the particle currents 
Ii ( t ) in a stationary situation, as quantified by the sym¬ 


metrized noise spectrum 

/ OO 

dTe^‘^^{SIi{T)6Ii,{0) + 6 Ii,{t)6Ii(0)).{25) 

-OO 

We use the MacDonald formula [lOl 01] together with the 
regression theorem [42] to obtain the two-time correlation 
functions {Sni{t)Sng{t + r)) for t —>• oo (cf. Appendix C) 
and find 


Sii 


H = 22 


1 +UJ' 


E 

k^O 


cos (k(Ri - R//)) 


r£/4 




(26) 


We have evaluated Eq. (26 1 for the infinite range interac¬ 
tion model, Eq. (20), and for the nearest-neighbor inter¬ 


action model Eq. (22) in d = 1 dimension. In both cases. 


the diagonal noise Z = Z' is independent of the channel 
index Z and displays scaling , i.e. a dependence on the 
dimensionless variable D = uj/{g'y) that in both cases 
contains the feedback strength parameter g, 


27 

5d=i 

27 


1 ^ A^-1 rz2 

N N 1 + rz 2 




2 I 


1 

1 -b 16/D2 


(27) 


(28) 


The noise spectra in these two cases have the correct 
high-frequency (or zero feedback g = 0) Poissonian limit 
S'(n —>■ oo) = 27 , cf. Fig. 

The opposite limit of zero frequency w = 0 corresponds 
to the Poissonian long-time Z —>■ 00 limit Eq. (21) at 


any finite N for Sod, whereas the N = 00 lattice model 
(nearest neighbor interaction) is determined by diffusive 
dynamics in d dimensions at small frequencies (that fails 
on short time- or large frequency scales). 


Sdjuj) _ 1 /^ 

27 2^2 V57 


0 , 


(29) 


which also follows from Eq. (23) with Eq. (26). The 


noise reduction at low frequencies highlights the strong 
suppression of fluctuations in each transport channel due 
to the feedback from all the other channels. 


IV. HEAT AND INFORMATION 
A. Condensation picture 

The interaction-induced feedback mechanism among 
the channels leads to a ‘condensation’ of the system into 
a state p(n, t) where fluctuations of the particle numbers 
ni are essentially suppressed at large times. 

Without feedback, the addition of the N indepen¬ 
dent Poissonian processes (channels) in the total parti¬ 
cle number ntot = jf (normalized by N) leads to 
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{Sntot{ty) = 'yt/N: the fluctuations of the macroscopic 
variable ritot are reduced by a factor 1 /N as compared to 
the fluctuations {5ni(tY) = in the individual, ‘micro¬ 
scopic’ channel variables ni. 

In contrast, with feedback the microscopic fluctuations 
are reduced at large times t, {Sni{t)‘^) —>■ jt/N, whereas 
the form Eq. ( |l^ guarantees that the fluctuations of ritot 
are not affected by the feedback: summation of Eq. (18) 
over I, use of = N6\^,q and the limit 

Ek^o = 0 again leads to (dntot(i)^) = jt/N, and the 
microscopic and macroscopic fluctuations are now of the 
same order. 

This explains the scheme shown in Fig. (Ic), where 
the channel variables ni are (continuous) position coor¬ 
dinates xi of a large cluster of N interacting Brownian 
particles on a line, with the cluster center-of-mass coor¬ 
dinate corresponding to ntot(t). The feedback (interac¬ 
tions) then condenses the (at time t = 0) ‘gas’ of initially 
independent particles into a tightly bound cluster, where 
fluctuations of relative distances between particles do not 
grow with time (as without feedback) but approach the 
fixed value obtained from Eq. as 

„ ,2\ 27^ l-cos(k(R,-Rp)) 

{{ni - nv) )t^oo = ^ 2^-E-■ (30) 


MO 


Tk 


B. Heat 

Starting from initially empty drains, this condensation 
process is accompagnied by an increase in interaction en¬ 
ergy £;(n), Eq. which has to be compensated by a 
flow of heat from the reservoirs that keep the drains at 
constant temperature and chemical potential all the time. 
Explicitly, a simple calculation leads to 




—Tit-t 


2/3 


(31) 


with the rates Tk defined in Eq. (19). By integration of 
Eq. (311, the quadratic interaction potential Eq. (|^ thus 
leads to the equipartition {E{n))oo — {E{n))o = ^ksTN, 
i. e. an interaction energy ^ksT stored in each channel. 

At the same time, the constant current I = 7 +— 7 _ > 0 
of particles (cf. Eq. @) generates a dissipative power 
—lAyi = I{yis — /id) > 0 in each channel that has to be 
dissipated into the reservoirs, and the total amount of 
heat flow Q into the reservoirs thus is 


Q = -NIAyi--{Ein))t. 


(32) 


We note that this simple splitting of Q into two contribu¬ 
tions is a consequence of our model for the rates Eq. (10), 


with its separation of excitation energies e;(n) = A/i; -|- 
AEi{n) into a constant single particle part A/i; and the 
part AEi{n) with the interaction energies Eq. (§. 

In the following, we check Eq. (321 by a somewhat more 


C. Fokker—Planck equation 

We use an approximation to the original rate equations 
Eq. where the integers ni become continuous, real 
number xi of which only Gaussian fluctuations are kept. 
In this way, we obtain a Fokker-Planck equation for the 
probability distribution P(x, t) in the usual form of a 
continuity equation 


N 


dr., ^ ^ d r, ^ 

-p(x,t) = - 


l=l 


dxi 


(33) 


7 d 


Ji{-K,t) = Fi(x)P(x,t) -P(x,t), (34) 

2 oxi 

with diffusion constant ^ and force term P/(x) = 7 _|_ — 
7 - +J2^=i9ij(.Xj-xi). 

Since Pi(x) is linear, this can be solved exactly and 
thus be used to derive explicit expressions for the Shan¬ 
non entropy 

S{t) =—{\nP)t = — J dxP(x, t) lnP(x, t). (35) 


Moreover, Eq. (33) is very convenient for directly reading 
the temporal chang 


off the splitting of the temporal change, 

d 
dt 


(36) 


into a (positive) entropy production rate n(t) and the 
global entropy flow 


N 




(37) 


from the total system into the reservoirs [13]. Crucially, 
even though S{t) and n(t) are ‘abstract’ quantities in the 
sense that they are no direct experimental observables, 
the entropy flow $(t) directly yields the (measurable) 
heat flow Q = 4d>(t) into the reservoirs. The expec¬ 


tation value Eq. (37) can be expressed in terms of the 
second cumulants J^iSxir) via the identification n; o xi 
from the Cw{t), Eq. , or from the direct solution of 
Eq. (33) (see below). Using Fourier transformation and 


some straightforward algebra, one obtains 

Q ^ l<i>(t) = tanh ^ - j^{E{n))t. (38) 

Now recalling that in the derivation of our classical tran¬ 
sitions rates after Eq. ^ we demanded /3|A/x| <gc 1, we 
may replace the tanh in Eq. (38) by its argument and 


elaborate argument based on entropy. 


thus recover our previous result Eq. (32) that was based 
upon a simple energy argument. 


D. Feedback information flow 

Finally, we quantify the information flow in the 
‘feedback-freezing’ of the distribution P(x, t). Eq. (33) 
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models the continuous counting variables xi as the posi¬ 
tions of N linearly coupled Brownian particles, a situa¬ 
tion analysed in detail (for iV = 2) by Allahverdyan et 
al. |44j via the mutual information T. Recently, Horowitz 
[5S] has generalized this to multipartite systems governed 
by Fokker-Planck equations like Eq. (33) and introduced 


an information flow based on the concept of neighbours 
influencing each other. 

Here, we use the simplest bipartite splitting of the 
whole system into a single channel I = 1 (with probability 
density Pi(xi, t)) and the remaining TV—1 channels (with 
probability density Pii{x 2 , where the mutual 

information 


I{t) = [ d^xP{x,t)\n — 

J Pi 


P(x,t) 




(39) 


quantifies the build-up of feedback effects within a small 
subsystem (single channel) embedded into a large feed¬ 
back environment. 

The temporal change i{t) = defines the 

two information flows between the two subsystems m 
133]. Using Eq. (33) and the equivalence of all channels. 


we can express the flows via entropies |44| . such as 

1 A 




(40) 


with the Shannon entropy S{t) of the total system, 
Eq. (35), and the quantity 


Si{t) = J d^x(^J^Ji{x,t)^lnPi{xi,t) (41) 

which has the form of a local Shannon entropy change, 
cf. Appendix C. 

In k-space, we find that the second cumulants crk(t) 
(defined in Eq. ( [To] )) completely determine the entropies 
via 

Slit) = ^4ln^crk(t), Sit) = i^^lncrk(t),(42) 


which in view of our Gaussian approximation is not so 
surprising, but demonstrates that there is a direct con¬ 


nection between the (abstract) information flow Eq. (40) 


and a fluctuation quantity that at least in principle is 
accessible via the counting statistics. This argument is 
in line with recent results by Ansari and Nazarov |45| . 
who derived a general connection between Renyi (and 
Shannon) entropy and full counting statistics. 

At small times t, owing to CTk(t —>■ 0) « jt, the in¬ 
formation flow is zero although both entropy changes 
are large, Siit) = jfSit) « the latter being due 
to the initial condition P(x, t = 0) = 6^ix). On the 
other hand, at any finite N and large times t, only the 
k = 0 component significantly contributes to Eq. (42) 
since CTk(t —>■ oo) « 5k, 07 ^ and the information flow de¬ 
cays towards zero again, iy^iit —>■ oo) ^ 

In its transient dynamics, thus, (.i'_>i(t) has a maxi¬ 
mum at some finite time of the order 517 , where g is the 


dimensionless coupling constant of the interaction model. 
This behavior is shown in Fig. ([^. 


For the infinite-range interaction model Eq. (20), the 


difference between local and total Shannon entropy van¬ 
ishes with increasing N since the notion of ‘boundary’ 
between subsystems becomes meaningless. Accordingly, 
with Slit) = JfSit) = gie~^^'^^ — 1)“^ for N ^ 00 , the 
information flow vanishes. On the other hand, at large 
but finite N, an interesting feature then is the plateau 
in 6 i'_>.i(t) at intermediate times, which is shown in the 
lower part of Fig. ([^ and which clearly corresponds to 
the ‘freezing’ of the fluctuations, i.e. the plateau in the 
second cumulants Cait). 


V. DISCUSSION 

The autonomous feedback control discussed here is 
based on a microscopic mechanism, but it comes with 


the constraints Eq. (11) and Eq. (12) for the interactions 
Uij, which suggests uncharged particles as appropriate 
candidates, such as ultracold atoms used in recent trans¬ 
port experiments |46| . In contrast, for charged particles 
(e.g., for electronic transport), the non-autonomous, syn¬ 
chronization feedback model Eq. (§ (which essentially is 
an open-loop control scheme in the space of particle num¬ 
bers n), appears to be more flexible. 

Needless to say that the choice of a particular control 
scheme depends on the specific control task, which in 
our model was to stabilize the full counting statistics at 
a given particle current I = 7 + — 7 _ over a longer time t. 
One aspect then is the issue of energy cost and efficiency. 
The two control schemes mentioned above both scale lin¬ 
early in t in terms of energy costs: the non-autonomous 
scheme by requiring a continuing modulation of energy 
barriers, and the autonomous scheme by needing large 
channel numbers N oc. t (and thus a large total heat flow 
jksTN oc t, cf. after Eq. (31)) if stabilization over large 
intervals t is required, as we saw in Fig. (2a). 

The analysis proposed here remained at the level of a 
Gaussian approximation for the full counting statistics, 
i.e. a linear Fokker-Planck equation. Our results might 
stimulate further work to explore the connection between 
the recently discussed flow of information in feedback net¬ 
works EH m na [3S], and the dynamical scaling type 
analysis for transport processes in interacting systems, 
such as in the Kardar-Parisi-Zhang (KPZ) model [35] for 
particle deposition and surface growth. For this latter 
point, the most urgent step then is to go beyond the lin¬ 
ear (in n) approximation Eq. (10) of the rates, either by 
going to higher orders in the expansion Eq. ([^ (i.e., be¬ 
yond the classical limit and towards lower temperatures), 
or by sticking to the classical limit Eq. 0 and employing 
interaction models beyond the simple quadratic model 
Eq. ([^. An interesting alternatively there could be the 
study of the KPZ equation (or other Langevin equations 
belonging to Eq. 0 ) per se from the perspective of en¬ 
tropy and information flow. 
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Appendix A: Non-autonomous feedback in 
single-channel model 


Here, we derive the equation of motion for the reduced 
probability 

pi{n,t) ='^Sn,,nP{n,t), (Al) 

n 

with n = (ni,nAf)^, for a single channel (chosen as 
I = 1 ) in the N ^ oo limit of the constant interaction 
model in Eq. ( [20| ). After multiplying Eq. ® with Sni,n 
and summing over n, only the I = 1 term oithe sum over 
I remains and 


■pi{n,t) (A2) 

= [7iT(n±ei)p(n±ei,t) -7i;p(n)p(n,t)] . 

n,± 


For gij = the rates have the form 7 i±(n) = 
fj. (ria^-v — rii ), with the average riav = (ni -I- ... -|- nN)/N 
and f±{x) = 7±(1 ± gx) for the linear rate model Eq. ([^ 
with feedback parameter g. Now, for N ^ oo the aver¬ 
age riav becomes a macroscopic variable. We neglect the 
fluctuations of the latter against those of ni by factoriz¬ 
ing 


(’^av - n) p{n, t) 

n n 

« h{{na-v)t-n)pi{n,t) (A3) 


and similar for the terms 7 i=p(n ± ei)p(n ± ei,t) in 
Eq. (A2). In this approximation, we effectively neglect 


the backaction of the first channel on all the other — 1 
channels. 

In this mean-field type approach, the effect of the other 
channels is that of an external ‘clock’ as expressed by the 
time-dependence of (nav)t = ( 7 + — 1-)^ = It according 
to Eq. (14). This leads to 




- f^{It-n)pi{n,t)], 


(A4) 


which describes non-autonomous feedback control as in 
[34] via rates Eq. ([^ with the ‘control current’ I = 
7 + — 7 _. As mentioned already after Eq. (21), the sec¬ 
ond cumulant Eq. ([^ of the model Eq. (A4) co-incides 
with that of the constant interaction iV-channel model, 
Eq. (21), for —>■ 00 . In that limit and within the 

Gaussian approximation, the correspondence therefore is 
exact. 


The symmetrized noise spectrum Sw^uj), Eq. (25), fol¬ 
lows from the MacDonald formula as outlined, e.g., in 
m- One starts from 


/ OO 

dre^^^Shit + T)5Iv{t) + (/ o /')), (Bl) 

-00 


which still depends on the initial time t that is sent to 
infinity at the end. We write 


j-t + T 

J dt'SIi{t') = Sni{t + t) - Sni{t), 


(B2) 


where Sni (t) = n; (t) — tli with tli = {ni (t)), cf. Eq. (14). 
As in im, this leads to 


Swjuj) 

2u} 


dTsin{uJT) 


(B3) 


^ ^ {{^'n‘i{t + T) - Sni{t)] [dni>{t + t) - 6ni>{t)]). 

To proceed, we now need (in the mixed terms in Eq. ( |B3| )) 
the two-time correlation functions {6ni{t)6nii{t + r)), 
which in view of Eq. (13) fulfill the regression equations 


dr 


{5ni{t)5ni,{t + r)) = 


(B4) 


N 


= ^^9Vj {{5ni(t)5nj{t + t)) - {j I')). 
1=1 


We can solve this by first transforming into Fourier (k) 
space similar to what we did to obtain the equal time 
correlation function Cw^t) = {5ni{t)Snii (t)) in Eq. ( |I^ , 
which in fact serves as the initial condition at r = 0 here. 
The result 


{5ni{t)6nii{t + t)) = (B5) 


N N 


E- 

k^O 


Fk 




can now be used in the MacDonald formula Eq. (B3), 
leading to 


nOO 

Siii{oj) = 2uj dTsin{ojT) 

Jo 


k^^O 


(B 6 ) 
■ (B7) 


Here, as usual [41] the sine function under the integral in 
the first term has to be understood as the lim it sinwr = 
lim„_>o and we now obtain Eq. (26). 
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Appendix C: Information flow 

From Eq. ( |39[ ), the temporal change of the mutual in¬ 
formation I{t) is 


N 


d 


i{t) = - E - In Pi - In Pp] , (Cl) 


1=1 


where we used Eq. ( p33| and the conservation of probabili¬ 
ties, J dxiPi{xi,t) = 0, f dX 2 ---dXNPl'{x 2 , = 0 , 

J d^xP(x,t) = 0. We can write 


d^ 


d ^ d \ 


dxi 

d 


1=2 


= / d^x ( ^ Ji(x,t) ) lnPi(a;i,t) 


(C2) 


by integration of parts of the second term, and similar 
for In Pi/. This yields the splitting j(f) = 6 i/_>i -I- ti_i.i' 
with the information flows 


ti'_j.i = — / d^x 


d 

dxi 


N 


Ji(x,t) In 


Pi 


P^i' = - J d^K\^^M^,t)yn^. (C3) 


Using the definitions of the global and local Shannon en¬ 
tropies, Eq. (35) and Eq. (41), we can now write the 


information flow ri'_i.i from channels 2, ..., N to channel 
1 in the form Eq. (40). 


For the evaluation of Eq. (41), we use integra¬ 

tion by parts to write 




(C4) 


We use the fact that Pi(xi,t) has to be a Gaussian, and 
we thus obtain (with the definition of Pi(x) in Eq. (33)) 
the explicit form 




N 




= 7 - 


El 


3 = 


2Ek^k 


= ^|lnE^k(t), (C5) 


where Sxi = xi — ( 7 + — 7 _)t and we used Fourier compo¬ 
nents, (Sxi) = ^ (Tk, and the second cumulants crk(t) 
in k space, Eq. ^1^ . 

For the global Shannon entropy change, in a similar 
manner (by exploiting the equivalence of all channels I in 
the homogeneous case) we obtain 


1 FkC 1 d 

k 

where again we used the fact that the full distribution 
function P is a Gaussian. 


2N ^ 1 

k 
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